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Abstract 

We present a quantum analysis of three-mode optoacoustic parametric interactions in an opti- 
cal cavity, in which two orthogonal transverse optical-cavity modes are coupled to one acoustic 
mode through radiation pressure. Due to the optimal frequency matching — the frequency sep- 
aration of two cavity modes is equal to the acoustic-mode frequency — the carrier and sideband 
fields simultaneously resonate and coherently build up. This mechanism significantly enhances 
the optoacoustic couplings in the quantum regime. It allows exploration of quantum behavior of 
optoacoustic interactions in small-scale table-top experiments. We show explicitly that given an 
experimentally achievable parameter, three-mode scheme can realize quantum ground-state cool- 
ing of milligram scale mechanical oscillators and create robust stationary tripartite optoacoustic 
quantum entanglements. 
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I. INTRODUCTION 



Optoacoustic interactions have recently become of great interest, for their potential in 
exploring the quantum behavior of macroscopic objects. Various experiments have demon- 

ed significantly through 



strated that the acoustic mode of a mechanical oscillator can be coo 
two-mode optoacoustic interactions 1Q|, [111, UJ: 
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14i | . The basic setup 



consists of a Fabry- Perot cavity with an end mirror. Linear oscillations of the mirror acoustic 
mode at frequency u> m scatter the optical-cavity mode (usually TEM o) at frequency ujq into 
Stokes (uq — u m ) and anti-Stokes (uq + uj m ) sideband modes which have the same spatial 
mode shape as the TEM o mode. The optical cavity is appropriately detuned such that the 
anti-Stokes sideband is close to resonance. Therefore, the anti-Stokes process is favored over 
the Stokes process. As a natural consequence of energy conservation, the thermal energy 
of the acoustic mode has to decrease in order to create high-energy anti-Stokes photons at 
uj + uj m . If the cavity- mode decay rate (related to the optical finesse) is smaller than the 
acoustic-mode frequency, theoretical analysis shows that these experiments can eventually 



achieve the quantum ground state of a macroscopic mechanical oscillator [15|, [16|, |17|, which 
would be a significant breakthrough in physics from both experimental and theoretical points 
of view. With the same scheme, many interesting issues have been raised in the literature, 
such as teleportation of a quantum state into mechanical degrees of freedom [3] , creation of 
stationary quantum entanglements between the cavity mode and the mechanical oscillator 
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20[ or even between two oscillators 2l|, |22|] . This in turn could be implemented in future 



quantum communications and computing. 

The concept of three-mode optoacoustic parametric interactions was first introduced and 
analyzed theoretically in the pioneering work of Braginsky et al. [l|. It was shown that 
three-mode interactions inside high-power optical cavities of large-scale laser interferometric 
gravitational-wave detectors have the potential to induce instabilities, which would severely 
undermine the operation of detectors. This analysis was elaborated by many other authors 
to more accurately simulate the real situation in next-generation advanced gravitational- 



wave detectors 
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26j and to find strategies for suppressing instability 
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281. 



Recently, the UWA group experimentally demonstrated three-mode interactions in an 80-m 
high-power optical cavity by exiting acoustic modes and observing resonant scattering of 
light into a transverse cavity mode 29]. 



2 



Different from the two-mode case, in three mode interactions, a single acoustic mode of the 
mechanical oscillator scatters the main cavity TEMoo mode into another transverse cavity 
mode which has a different spatial distribution from the TEM o mode. Specifically, when 
the TEM o mode is scattered by the acoustic mode, the frequency is split into Stokes and 
anti-Stokes sidebands at Uq ± u m , and in addition, the spatial wave front is also modulated 
by the acoustic mode. Three-mode interactions happen strongly when both the modulation 
frequency and spatial mode distribution are closely matched to those of another transverse 
optical-cavity mode. Under these circumstances, both the carrier and sideband modes are 
simultaneously resonant inside the cavity and get coherently build up. Taking into account 
the resonance of the acoustic mode, the system is triply resonant, with the interaction 
strength scaled by the product of the two optical quality factors and the acoustic quality 
factor. If the transverse optical cavity mode has a frequency lower than the TEMoo mode, 
the Stokes sideband will be on resonance and the interaction provides positive amplification 
of the acoustic mode, while if the transverse mode has a frequency above the main cavity 
mode TEMoo, the anti-Stokes sideband will resonate, and the system has negative gain and 
the acoustic mode will be cooled. The underlying principle of both two and three mode 
interactions is similar to the Brillouin scattering, except that the modulation occurs not 
through changes in refractive index of the medium, but through bulk surface motion of a 
macroscopic mechanical oscillator (i.e. the acoustic mode) which modulates the optical path 
of the light. 

While three-mode interactions are inconvenient byproducts of the design of advanced 
gravitational-wave detectors, they can be engineered to occur in small scale systems with 
low mass resonators, which can serve as an optoacoustic amplifier and be applied to acoustic- 



mode cooling 



30, 



32|. Besides, due to its triply resonant feature, the three- mode system has 



significant advantages compared with the two-mode system and allows much more stronger 
optoacoustic couplings. To motivate experimental realizations, we have suggested a small- 
scale table-top experiment with a milligram mechanical oscillator in a coupled cavity jsoj]. 
Using the extra degree of freedom of the coupled cavity, the cavity mode gap ( the difference 
between the two relevant cavity modes) can be continuously tuned such that it is equal to 
plus or minus the acoustic-mode frequency, which maximizes the three-mode interactions 
strength. We also pointed out that, in the negative-gain regime, this experimental setup 
can be applied to resolved-sideband cooling of a mechanical oscillator down to its quantum 
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ground state. In that paper, we used the classical analysis presented by Braginsky et al. to 
obtain the effective thermal occupation number n of the acoustic mode. This analysis breaks 
down when n < 1 and the quantum fluctuations of the cavity modes have to be taken into 
account. To overcome this limitation, we used the similarity in the Hamiltonian of the two- 
mode and the three-mode system, and argued that the quantum limit for cooling in both 
systems is the same without investigating the detailed dynamics. However, in order to gain a 
quantitative understanding of the three-mode system in the quantum regime, it is essential 
to develop a full quantum analysis which includes the dynamical effects of the quantum 
fluctuations. Besides, as we will show, the quantum analysis reveals a most interesting 
non-classical feature of three-mode systems: stationary tripartite quantum entanglement. 

The outline of this paper is the following: In Sec. [Ill we start from the classical analysis 
given by Braginsky et al. and then quantize it with the standard approach. In Sec. IHH 
we use this quantized Hamiltonian as the starting point to analyze the dynamics of the 
three-mode system. Further, based upon the Fluctuation-Dissipation-Theorem (FDT), we 
derive the quantum limit for the achievable thermal occupation number in cooling experi- 
ments. To motivate future small-scale experiments, we provide an experiment ally- achievable 
specification for the quantum ground-state cooling of a mechanical oscillator. In Sec. IIVI 
we investigate the stationary tripartite optoacoustic quantum entanglement and we show 
that the same specification for the cooling experiments can also be applied to realize robust 
stationary optoacoustic entanglements. 



II. QUANTIZATION OF THREE-MODE PARAMETRIC INTERACTIONS 



brmulations of three-mode optoacoustic 
l] and then apply standard procedures 



In this section, we will first present the classical 
parametric interactions given by Braginsky et al. 
to obtain the quantized version. 

Classical Picture. A detailed quantitative classical formulation of three-mode interactions 
was given in the Appendix of Ref. A Lagrangian formalism was used to derive the 

classical equations of motion and analyze the stability of the entire three-mode optoacoustic 
system. The formalism can be easily converted into Hamiltonian language, which can then 
be quantized straightforwardly. For convenience, we will use slightly different notation and 
definitions for the optical fields. Further, we assume that the two optical-cavity modes are 
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FIG. 1: Spatial shapes of the TEMoo an d TEMoi modes and the acoustic torsional mode. 

the TEMqo an d TEM i modes and the acoustic mode has a torsional mode shape (about 
vertical axis) which has a large spatial overlap with the TEMoi mode as shown in Fig. 
[TJ It can be easily extended to general cases with other transverse optical modes and 
acoustic modes. Assuming the electric field is linearly polarized in the transverse direction 
perpendicular to z axis, the electromagnetic fields (E, H) of the cavity modes can be written 
as 

EM = (^j f t (r ± )sm(k t zMt), (1) 
Hi(t) = T (^t) ' /i(rL)cos(M)ft(t). (2) 



h \e V / 

Here i — 0,1 represent the TEM o and TEM i modes; fi(r±) are the transverse mode shapes; 
Ui denote the eigenfrequency; ki are the wave numbers; V is the volume of the optical 
cavity; <&(t) are the generalized coordinates of the fields; qi(t) are the time derivatives of 
<7i(£). At present stage, the appearance of huji is just to make the generalized coordinates 
dimensionless. The classical Hamiltonian of this system is given by 

H = H m + ~J dr ± (L + xu z )[e {E + E t ) 2 + fi (H + #i) 2 ], (3) 

where L is the length of the cavity; x is the generalized coordinate of the acoustic mode; u z 
is the vertical displacement. The free Hamiltonian of the acoustic mode is 

H m = -hu m {<fi n + p} n ) (4) 

with q m = xj y/hj (muj m ) and p m is the momentum normalized with respect to \fhmu m . 

After integrating over the transverse direction and taking into account of the mode shapes, 
we obtain 

H = H m + H0 + H1 + H mt . (5) 
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Denning dimensionless canonical momentum Pi(t) = qi(t)/uji, the free Hamiltonian of the 
two cavity modes are 

'Hi= l -huj l {q 2 l +Pl) (6) 
and the interaction Hamiltonian is given by 

Hint = hGoqmtqoq! + p pi), (7) 

where the coupling constant is defined as Go = \/A h ujqU\ / (m u} m L 2 ) with the geometrical 
overlapping factor A = (L J dr^Uzfofi/V) 2 . 

Given above Hamiltonian, it is straightforward to derive the classical equations of motion 
and analyze the dynamics of the system, which would be identical to those in the Appendix 



of Rcf . 



1|. To quantify the strength of three- mode interactions, Braginsky et al. introduced 



the parametric gain 1Z, as defined by 

K = ± 2AI ^ = ± 2AIoQoQ 1 Q I n , g) 
mu; m L 2 7o7i7 m mujQUJ^L 2 

where ± correspond to either positive gain or negative gain; Jo is the input optical power of 
the TEMoo mode; we have defined optical and acoustic-mode quality factors Qi = u)i/^i{i = 
0, l,m). Due to optoacoustic interaction, the decay rate 7 m of the acoustic mode will be 
modified to an effective one j' m , which is 

-rfXl-tt) 7ra . (9) 

When 1Z > 1, the decay rate becomes negative and this corresponds to instability. In this 
paper, we are particularly interested in the regime where TZ < which gives rise to the 
acoustic-mode cooling. The effective thermal occupation number n[ h of the acoustic mode 
is given by 

_/ ^th7m n t h . . 

= — = T^i (10) 

with nth denoting the original thermal occupation number. It looks as if — 1Z ^> 1, n' th can be 
arbitrarily small. However, in this case, the classical analysis breaks down and the quantum 
fluctuations of the cavity modes will set forth a quantum limit for the minimally achievable 
n' th , which will be detailed in the following quantum analysis. 

Quantized Hamiltonian. The classical Hamiltonian derived above can be quantized by iden- 
tifying these generalized coordinate and momentum as Heisenberg operators which satisfy 
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following commutation relations 

[qj,Pj>]=i6jf, (j,f = 0, l,m). (11) 
The quantized Hamiltonian is then given by 

H = ^ + Pi) + hGogmiMl + PojPl) + W e xt, (12) 

i=m,0,l 

where we have added 7i cxt to take account of the coupling between cavity modes and external 
continuum optical fields due to the finite transmission of the cavity. This Hamiltonian is 
convenient for discussing stationary tripartite quantum entanglement as will be shown in 
Sec. |IVj for those generalized coordinates and Pi correspond to the amplitude and phase 
quadratures in the quantum optics entanglement experiments. 

To discuss the ground-state cooling as will be investigated in Sec. IHIt it would be 
illuminating to introduce annihilation operators for the two cavity modes a = (qo + ipo)/y/2 
and b = (q± + ipi)/y/2, such that the normally ordered quantized Hamiltonian can be 
rewritten as 

H = -^w m {f m + p 2 m ) + hw a^a + HuiPb + hG Q q m (a ] b + b ] a) + H cxt . (13) 



III. QUANTUM LIMIT FOR THREE-MODE COOLING 

In this section, we will start from the Hamiltonian in Eq. (1131) to derive the dynamics 
and discuss the quantum limit for the ground-state cooling experiments using three-mode 
optoacoustic interactions. As we will see, due to similar mathematical structure as in the 
two-mode case, the corresponding quantum limit for three-mode cooling is identical to the 



15| and Wilson-Rae et al. [3] in the 



resolved-sideband limit derived by Marquardt et al. 
two-mode case. 

Equations of Motion. The dynamics of this three-mode system can be derived from the 
quantum Langevin equations (QLEs). In the experiments, the TEMoo mode is driven on 
resonance at ujq. Therefore, we choose a rotating frame at uq, obtaining the corresponding 
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nonlinear QLEs as: 

q m = UJ m p m , (14) 

p m = -uj m q m - 7 m p m - G (a ] b + Pa) + f th , (15) 
a = --f a-iG q m b+ ^/2rfyd in , (16) 
b = _( 7l _ i A)6 - iG q m a + ^/2r/lb in . (17) 

Here the TEM o and TEM i mode gap is given by A = oj\ — ujq; Go(a>b + ¥a) corresponds to 
the radiation pressure which modify the dynamics of the acoustic mode and is also respon- 
sible for the quantum limit; we have added thermal noise £ t h whose correlation function, 
in the Markovian approximation, is given by (£th(£)£th(0) = 27 m n t h5(t — t'). In obtaining 
above equations, we have also used Markovian approximation for 7i ex t, namely 

n cxt = % h{^2^a)a m + v^n 6 f 6in - H.c.) (18) 

with H.c. denoting Hermitian conjugate. 

To solve the above equations, we can linearize them by replacing every Heisenberg oper- 
ator with the sum of a steady part and a small perturbed part, namely 6 = o + Sd(e) with 
e<l, We treat Brownian thermal noise £ t h, the vacuum fluctuations ^/yoSd in , y/yi.5b- m and 
Sq m as the order of e. In the experiments, the TEM o mode is pumped externally with a 
large classical amplitude a- m while the TEM i mode is not with b m = 0. Therefore, to the 
zeroth order of e, the steady part of the cavity modes are simply given by 



a = ^/VToOin = vW(7oMi b = -iG aq m . (19) 

Without loss of generality, we can set q m = 0. Therefore, b = and this allows us to 
eliminate the TEM o mode from the first-order equations, which are 

5q m = u m 5p m , (20) 

Spm = -u m 8q m - 7m $Pm - G a(5b + 5&) + £ th , (21) 

6b = -(71 + i A)5b -iG a5q m + y/^6b- m . (22) 

Here we have chosen an appropriate phase reference such that d- m is real and positive. The 
above equations can be solved in the frequency domain, namely 

oU£p(n) + &*(»)] 

qm[il) ~ ~(n*-u&) + i7n.fi' [ } 

^/^n = G g ^gm(fi) + Zy / 27i'^in(fi) ^ 

(fi - A) + 271 
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^~ Wm+ f.. a\9 i _.2ir/., , A^i^T ^'J 



where the radiation pressure 

~ _ 2G 2 a 2 A 5g m (ty ~ 2G aV7r[( 7 i - in)S^(Q) - A gp 2 (fl)] 

rpl >~~ [(n-A) + i 7 i][(n + A) + <7i] 1 j 

with amplitude and phase quadratures Sq2 (fi) = [<J&(fi)+<Sir(— fi)]/V2 and 5f>2(^) = [5&(^) — 
f2)]/(\/2z). In the expression of F Tp , the part proportional to Sq m is the optical spring 
effect. For a high quality-factor oscillator with u m j m , the decay rate 7 m and the 
eigenfrequency u m of the acoustic mode will be modified to new effective 7^ and u' m , as 
given by 

y =7 + 4Gga 2 A^ m7l 

7m 7m [(« m - A)^ + 7 2][( Wm + A)^ + 7?] ' ^ } 

G 2 a 2 A(u 2 m - A 2 - 7l 2 
[(^ m -A) 2 + 7 1 2 ][(^ m + A) 2 + 7i 2 
In our case, the TEM o and TEM i mode gap is A = uj\ — u = u m . For the resolved- 
sideband with 7x <C u m , we obtain 

G 2 a 2 G 2 a 2 
7™~7 m + — ; w^«w ro --^-. (28) 

7l 4 ^m 

If we define the parametric gain as 71 = (j m — 7 / t )/7 m , then in this case 

n _ G la 2 _ 2A/ o^i (29) 
7i7m ma; m L 2 7 7i7 m ' 

This is identical to Eq. (JSj) in the negative-gain regime, which was obtained from classical 

analysis by Braginsky et al. However, different from Eq. (fTUj) . the resulting thermal 

occupation number of the acoustic mode is given by 

-/ ^th7m . _ n th 

n th = — ; h n quant = — + n quant (30) 

7m I -71 

where the extra term n quant originates from the vacuum fluctuations in F Tp , i.e. terms 
proportional to 5p2 and 5q~2- Since in the case of large 7Z or equivalently strong optoacoustic 
coupling, n' th pa n quant and the acoustic mode will finally reach a thermal equilibrium with the 
cavity modes. The lowest achievable thermal occupation number n quan t will be determined 
by this optical heat bath (cavity mode + external continuum mode). 

To derive this quantum limit n qU ant, we will apply the Fluctuation-Dissipation-Theorem 
(FDT). Specifically, given any two quantities A(t) and B{t) which linearly depend on field 
strength, we can define the forward correlation function 

C M (t - if) = (A(t)B(t>)) (t>t>), (31) 



where ( ) denotes the ensemble average. According to the FDT, we have 

S M (Q) + S M (-Q) _eP hn + l 



or equivalently, 



2n eff (0) + 1, (32) 



Mn) = sj^s^ry (33) 

where 5^(0) is the power spectral density (Fourier transform of C^b) an d P = l/(^B^eff) 
and the effective thermal occupation number n c g = 1/ (e l3hn — 1). In our case, we can simply 
substitute A, B with the amplitude of the TEM i mode Sb by fixing q m = 0. From Eq. (|24p 
and using the fact that for vacuum fluctuation (56^(0)56^(0')) = 27r5(0 — O'), we obtain 

■%um = (n _^, +7 , - (34) 

Since the acoustic mode have a very high intrinsic quality factor (cu m ^> j m ), the energy 
transfer between the cavity modes and the acoustic mode only happens around u m . There- 
fore, from Eq. (1331) and Eq. (1341) . the final quantum limit is given by 

X 

where we have used the fact that for the resonant case, A = u)\ — loq = u m . To achieve the 
quantum ground state, i.e. n quant ~ 0, we require u m 3> 71 and this is simply the resolved- 



riquant ~ n cS (u m ) = [ ) , (35) 



sideband limit obtained in the pioneering works of Marquardt et al. 15|] and Wilson-Rae et 
al. 3- 

The reason why the quantum limit for three-mode cooling is identical to the two-mode 
case can be readily understood from the fact that the TEM o mode is eliminated from the 
optoacoustic dynamics as shown explicitly in Eq. (12171) (T2~T|) (T2"2"|) and we essentially obtain an 
effective two-mode system. As suggested by Chen 31|, this equivalence can be more obvious 



by mapping this three-mode system into a power and signal-recycled laser interferometer, 
as shown in Fig. [2j The TEM o and TEM i modes can be viewed as the common and 
differential modes in the interferometer respectively. The torsional mode corresponds to 
the differential motion of the end mirrors and A is equivalent to the detuning of the signal 
recycling cavity. In the power and signal-recycled interferometer, even though there is 
no high-order transverse optical mode involved, the two degrees of freedom of the power 
recycling mirror and the signal recycling mirror enable simultaneous resonances of the carrier 
and sideband modes, which is achieved naturally with the three-mode optoacoustic scheme. 

10 




FIG. 2: Equivalent mapping from three-mode system to an power and signal-recycled interferom- 
eter. The TEMoo and TEMoi mode can be viewed as the common and differential optical modes 
in the interferometer respectively. The torsional acoustic mode is equivalent to the differential 
motion of two end mirrors in the interferometer. By adjusting the positions of power-recycling 
mirror (PRM) and signal-recycling mirror (SRM), we can make the carrier and sideband modes 
simultaneously resonate inside the cavity, the same as the three-mode scheme. 



The above discussion shows that mathematically two-mode interactions and three-mode 
interactions are very similar. However, it is very important to emphasis that from an ex- 
perimental point of view, there is an important difference. Specifically, the steady-state 
amplitude a in the radiation pressure F rp is amplified by the optical resonance while for the 
two-mode case this amplitude is highly suppressed due to large detuning. In other words, in 
order to achieve the same optoacoustic coupling strength experimentally, the input optical 
power in the two-mode scheme needs to be 1 + (A/70) 2 times larger than the three-mode 
scheme. This is a large factor in the resolved-sideband regime with A ^> 70 (the opti- 



mal A = u m ). Besides, in three-mode interactions, the condition A = uj m e 
optimizes the energy transfer from the acoustic mode to the cavity mode 
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so naturally 



171 ]. There- 



fore, three-mode scheme greatly enhances the optoacoustic coupling and is able to achieve 
resolved-sideband limit without compromising the intra-cavity optical power. As mentioned 
in Ref. 32|, the amplitude and laser phase noise can also be reduced significantly with 



three-mode scheme, simply due to the filtering of the cavity resonance. 

To motivate future cooling experiments with three-mode interaction, now we present an 
experimentally achievable specification for the quantum ground state cooling of a milligram 
scale mechanical oscillator. We choose that the mass of the mechanical oscillator m = 0.1 mg; 
the length of the cavity L = 2 cm; the acoustic-mode frequency u m /2-ir = 10 6 Hz, the 
acoustic-mode quality factor Q m = u rn /'~f rn = 10 7 ; the optical finesse T = 10 4 . Given an 
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input optical power of the TEMoo mode Jo = 50 mW and the environmental temperature 
T = 4 K, the corresponding effective thermal occupation number of the mechanical oscillator 
~ 0.5. 



IV. STATIONARY TRIPARTITE OPTOACOUSTIC QUANTUM ENTANGLE- 
MENT 



As shown in the works of Vitali et al. [l9( and Paternostro et al. 20], optoacoustic 
interaction provides an very efficient way of generating stationary quantum entanglements 
among cavity modes and the acoustic mode. Once experimentally realized, it will have 
significant impacts on future quantum communications. Following their formulism, we will 
investigate the stationary tripartite quantum entanglement in the three-mode optoacoustic 
system by first analyzing the dynamics and then evaluating the entanglement measure - 



logarithmic negativity E^- defined in Ref . [33|, [34 ] . 



Starting from Hamiltonian in Eq. (|T2|) . the corresponding nonlinear QLEs in the rotating 
frame at the laser frequency lol can be written as 

Qm = ^mPm, (36) 

p m = -uj rn q m - 7 m p m - G (<?o<?i + PoPi) + £th, (37) 

% = -7o9o + AoPo + G^qmpx + v^7b9o n > ( 38 ) 

Po = -lo'Po ~ A g - G q m qi + V^ToPo 1 ' ( 39 ) 

O-i = -7i<?i + Aipi + G q m p + v^gl 11 , (40) 

Pi = -7iPi - Aigi - G q m q + (41) 

where A = ujq — uj^ and Ai = lo\ — uj^. Slightly different from the cooling experiments, 
here we need to externally drive both the TEMoo and TEMoi mode simultaneously to create 
tripartite quantum entanglement. We choose an appropriate phase reference such that the 
classical amplitude pi = and qi ^ (z = 0, 1) which is related to the input optical power 



Ii by qi = \f2Ii/ (hw^i). Similar to the previous case, we can linearize above equations as 



x 



r -Mx T + n T , (42) 
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with T denoting a transpose transformation and 

= ( 8q m , 5p rn , Sq , 5p , Sq u 5pi 
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( o, Cth, 

and matrix M is given by 



(43) 
(44) 



M 









V G % 



\ 

7 m G q 1 G q 

-7o A 

-A -7o 
-7i Ai 
-Ai - 7l J 



(45) 



At first sight, the mathematical structure is identical to the one analyzed by Paternostro 



et al. 



20|. Apart from differing in the coupling constants (here we need to consider the 



overlapping factor A), there is another important difference: After linearization, the ra- 
diation pressure term G Q (%qi + PoP i) in Eq. (j3~7|) is proportional to %5q\ + q\5% rather 
than q q — q^i considered in Ref. 20(. As we will show, similar to the case for cooling 
experiments, the coherent build up of both the TEM o and TEM i mode and optimal mode 
gap <jj\ — c^o = enhance the entanglement significantly, which make it easier to achieve 
experimentally. 

Assuming the system is stable, i.e. all eigenvalues of M have negative real part, the 
stationary solutions to Eq. (1421) can be written down formally as 



POO 

^(oo) = V / dt'ie^'^h^t') 

A JO 



(46) 



where we have neglected the initial-condition terms which decay away as the system ap- 
proaches the stationary state. We assume that all the noises are Markovian Gaussian pro- 
cesses and the correlation functions are 



(Tijit-t 1 ) = DySit-t') 



(47) 



where are the elements of matrix D and = Diag[0, 2j m kBT/(?kd rn ), 7o,7o>7i) 7i]- 
The corresponding stationary covariance matrix among the cavity modes and the acoustic 
mode can then be written as 



POO 

V(oo) = / dt[e Mt ]B[e Mt ] T , 
Jo 



(48) 
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and the components of V can be obtained by solving following algebra equations: 



MV + VM 1 = -D. 



(49) 



For this tripartite continuous- variable system (one acoustic mode + two cavity modes), 
one necessary and sufficient condition for separability is the positivity of partially transposed 
covariance matrix 35, 0, [3?]]. In our case, partial transpose is equivalent to time reversal 
and can be realized by reverting the momentum of the acoustic mode form p m to — p m , 
namely 



V 



pt 



VI 



(50) 



By evaluating the positivity of the eigenvalue of V pt , we can directly determine whether 
entanglement exists or not. To reveal the richness of the entanglement structure, we will 
not directly analyze the positivity of V„ for the entire system, but rather following Ref. 
[201 ]. we look at the entanglement between any bipartite subsystem using the logarithmic 
negativity E^. Given the 4x4 covariance matrix V su b for any bipartite subsystem, 



V s 



sub 



^2x2 ^2x2 
^2x2 B 2X 2 



(51) 



the logarithmic negativity E^f is defined by 34j, [35] 



Ej^j = max[0, — In 2<r_ 



(52) 



with a- = - ^/W^- 4 det W mh /y/2 and £ = det A + det B - 2 det C. 

For numerical estimations, we will use the same specification as given in the previous 
section for the cooling experiments. We will focus on the situation relevant to the experi- 
ments with oj\ — ujq = uo m and the TEMoo mode driven on resonance (Ao = 0, Ai = uo m ). In 
Fig. [3j we show the resulting E_\f as a function of the input optical powers of both optical 
modes. Given the specifications, the entanglement strength between each optical mode and 
the acoustic mode becomes stronger as the optical power of their counterpart increases (un- 
til the system becomes unstable). This is understandable, because we have q m (qoqi + qoQi) 
type of interaction, and the coupling strength between the TEM o mode and the acoustic 
mode directly depends on the classical amplitude of the TEM01 and vice versa. For the en- 
tanglement between two optical modes, it reaches maximum when both optical modes have 
medium power. This can be attributable to the fact that the entanglement between the two 
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FIG. 3: Logrithmic negativity Ej\f as a function of the input optical powers of both modes. Other 
specifications are identical to those for cooling experiments given in the previous section. The 
left panel shows Effi for the entanglement between the TEMoo mode and the acoustic mode; The 
middle panel presents Effi for the TEMoi mode and the acoustic mode; The right panel shows 
E$ for the TEM 00 mode and the TEM i mode. 

optical mode is mediated by the acoustic mode, and both Effi and Ej^ 1 should be large to 
give a reasonable E^. Besides, as shown explicitly in Fig. HI the condition to\ — u = oj m 
will naturally optimizes the entanglement between the TEM i mode and the acoustic mode. 
This is because the Lorentzian profiles of the TEMoi mode and the acoustic have the largest 
overlap when A = u m . In this case, both the TEMoi mode and the acoustic mode are 
driven by the same vacuum field, which gives the maximal entanglement. Therefore, the 
optimal condition for the cooling experiment will simultaneously optimize the entanglement 
strength, as also been observed by Genes et al. [38]. 

To illustrate the robustness of this tripartite entanglement, we show the dependence of 
Ej^f on the environmental temperature in Fig. [5j The entanglement between the optical 
modes and the acoustic mode is very robust and it persists even when the temperature goes 
up to 80 K. Although the entanglement between the two optical modes is relatively weak, 
yet it changes slower as the temperature increases, and it vanishes when the temperature 
becomes higher than 15 K. The robustness of the optoacoustic entanglement was also shown 
previously by Vitali et al. [3]. This is attributable to the strong optoacoustic coupling 
which suppresses the thermal decoherence of the acoustic mode. With both the TEMoo 
mode and the TEMoi mode on resonance, we can obtain much higher intra-cavity power 
compared with the equivalent detuned two-mode system. Given moderate input optical 
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FIG. 4: Logarithmic negativity Ej^ 1 as a function of cavity modes gap A = u\ — ujq. As we can 
see, the condition A ~ u m , which optimizes the cooling, also maximizes the entanglement between 
the TEMoi mode and the acoustic mode. Here we have assumed Iq = 4.5 W (Higher Iq will make 
the system unstable for small A) and I\ = 0. W. Since it can be viewed as an effective two-mode 

n 

system in this case with I\ = 0, it simply recovers the results given by Vitali et al. [191 ]. 




FIG. 5: Logarithmic negativity as a function of temperature. The solid curve stands for Effi, the 
dashed curve for Effi and dash-dot curve for Ejj . We have chosen the optimal parameters for each 
curve. 



power, this allows us to achieve stronger entanglement between the optical modes and the 
acoustic mode of a massive mechanical oscillator (~ mg). Of course, this robustness of 
entanglement is conditional on the fact that the mirrors of the cavity can sustain a high 
optical power ~ 10 4 W. If the beam size is of the order of mm, this corresponds to a power 
density of around 10 6 W/cm 2 , which is achievable with the present technology 39]. 

To verify this tripartite entanglement experimentally, we can apply the same protocol as 



proposed in Ref. 
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Specifically, through measuring the outgoing field, we can 



build up statistics and construct the covariance matrix V cxp of this tripartite system based 



1(3 



on the measurement results and then analyze whether the partially transposed covariance 
matrix VP* p fails to be positive definite. If VP* p has a negative eigenvalue, this will give an 
unambiguous signature for quantum entanglement, because any classical correlation always 
gives a positive definite V^* . Besides, we can also use V cxp to evaluate the logarithmic 
negativity of any bipartite subsystem to determine whether entanglement exists or not 
in a given subsystem. Since the tripartite entanglement is stationary, this means that the 
optoacoustic interactions protect the quantum entanglement from the thermal decoherence 
which is a significant issue in non- stationary quantum entanglements. In principle, we can 
make a sufficiently long integration of the output signal such that the shot noise is negligibly 
small and V exp would be a direct verification of what we have obtained theoretically. 

V. CONCLUSION 

We have analyzed the three-mode optoacoustic parametric interactions in the quantum 
picture. We have derived the quantum limit for the cooling experiments with three-mode 
interactions based upon the Fluctuation-Dissipation-Theorem. Besides, we have shown the 
existence of the tripartite quantum entanglements in this system. The simultaneous reso- 
nances of the carrier and sideband modes in the three-mode system allows more efficient 
acoustic-mode cooling and more robust optoacoustic entanglement than the two-mode sys- 
tem. This work provides the theoretical basis for the feasibility of realizing both ground-state 
cooling and stationary optoacoustic quantum entanglements using three-mode optoacoustic 
parametric interactions in small-scale table top experiments. 
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